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Radar offers promise of orbifal monitoring of near-surface soil moisture* The 
complex dielectric constant affords an excellent estimate of the soil moisture content 
due to its dependence of the free water content in the soil , Target composition 
dictates the complex dielectric constant while surface roughness is characteristic 
of surface and to some degree sub-surface composition. By varying any one of the 
controllable parameters of a radar system (incidence angle, polarization, frequency) 
it is possible to measure the effects of surface roughness and dielectric constant from 
the standpoint of returned power. The effect of roughness on the radar backscatter- 
ing coefficient, can be minimized by the proper choice of the radar parameters 
without any significant reduction in the sensitivity to variations in soil moisture. 

Soil moisture is a greatly fluctuating entity and therefore monitoring must be per- 
formed relatively often and with realistic sensor instrumentation to insure optimum 
data retrieval. The following recommendations are submitted: 


angle of incidence range; 7° to 15° 

frequency: 4 GHz 

polarization: HH or W 

revisit time; 4 days 

resolution: 100 meters 


• « • 
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RADAR AND RADIOMETER MEASUREMENT OF SOIL MOISTURE — 

STATE OF THE ART 

Brad Hanson 

1.0 INTRODUCTION 

The surface layer of soil is of utmosf importance to civilization. Maximizing 
cultivation output relies upon the maintenance of an optimum environment for this 
layer, o goal of cultivators for centuries. A key parameter in this environment is 
water but it has been only recently that changes in the soil water content for this 
surface zone have been understood [Jackson, 1973], 

The purpose of this paper is to review recent work in the field of remote 
sensing relative to soil moisture, to recognize the target parameters necessary if 
optimum data retrieval is to be realized, and to recommend proper sensor instru- 
mentation to achieve this goal. 

2.0 PREVIOUS WORK 


2,1 Passive Sensors 

The primaiy concern to the hydrologist Involved in large-scale, water resource 
management in farming regions is soil moisture; the method most actively utilized for 
data gathering lies in the realm of remote sensing, A number of investigations to deter 
mine soil moisture have been undertaken to determine the potential use of remote 
sensors which operate in the optical and thermal infrared regions of the electromag- 
netic spectrum [Luder, 1959; Gates, 1964; Stockhoff, 1971; Werner et al., 1971; 
and Reeves, 1973], Their results reflect two major problems; cloud cover often 
hampered data retrieval and the sensor response at optical and infrared frequencies 
produced information relative to only a very thin layer at the air-soil interface. 

Extensive research with passive microwave radiometers has been undertaken 
by many institutions including pr' .'ute industry, government agencies, and various 
universities. 


Edgertori [1968] conducted experiments utilizing a microwave radiometer to 
ascertain its app!lcc6ility for determining soil moisture. He noted that differences 
in emission could be measured from fields exhibiting variations in soil moisture con“ 
ditions. Similar results were noted by Poe [19711. Schmugge and others [19721 
concluded that it was possible to monitor soil moisture from a radiometer mounted 
on an aircraft at the frequencies of 1.42 GHz, 4,99 GHz, 19.35 GHz, and 37 GHz. 
They noted, however, that the emissions were not only a function of soil moisture, 
but were affected by soil type and surface conditions such as roughness and vegeto“ 
tive cover. Experimentation by other research groups [Jean et al., 1972; Siu, 

1974] substantiate the influence thot soil moisture possesses on radiometric emissions. 

There exist two major drawbacks which place severe limitations on an op era “ 
tional system for continuously monitoring soil moisture. Atmospheric effects such 
as cloud cover would severely limit or negate data, thus placing restrictions in 
some instances on the usefulness of the information retrieved. Achievable resolution 
is another problem. Edgerton [1968] and Poe [19711 noted that at satellite altitudes, 
soil moisture discrimination is practical only for gross spaciai differences; at altitudes 
feasible from aircraft, the swath width involved is much too narrow, A similor 
problem Was noted for a radiometric study over a Texas site [Sobti, 1975] . He 
concluded that while soil moisture has a substantial influence on the radiometric 
response, it is extremely difficult to predict soil moisture using the gross resolution 
sensors. It was noted, however, that it was feasible in a few instances to refine 
estimates where crude ones existed. 

Operational constraints, both from the standpoint of ^stem parameters and 
target response, severely limit the amount of useable data possible from passive 
microwave radiometers. Resolution is a major concern along with atmospheric 
restrictions such as cloud cover. It is the opinion of this investigator that launch- 
ing an operational satellite system to monitor soil moisture with a passive microwave 
radiometer as the principle sensor would yield a limited amount of useabie data. 
However, a satellite incorporating an active microwave sensor such as a radar plus 
o radiometer may lead to more information about a given target than if the latter 
system were eliminated from a satellite, sensor package. 

■ 2.2 Active Sensors 

Active remote sensors such qs radar are not restricted by the problems pre" 


viously discussed. Radar is nof hampered by non-precipifafion clouds and depending 
upon frequency, surface as well as near“sub-surface moisture condition information 
is practical. Also, fine resolution imagery is achievable from any altitude since 
resolution is independent of range (fully focussed, synthetic-aperture radar). 

Five parameters affect the reflection or scattering of electromagnetic waves 
from a target; they are polarisation, frequency, incidence angle, complex dielectric 
constant and surface roughness. Of these parameters, the last two are controlled by 
the target. 

Radar backscatter is dependent upon not only the surface and sub-surface 
geometry, but also upon the dielectric properties. The energy incident on the 
terrain surface is "specularly" or "diffusely" reflected in varying proportions de- 
pending upoi3 the roughness of the terrain, Surface roughness is a geometric prop- 
erty of the terrain; it is not an absolute roughness, but rather roughness expressed 
relative to wavelength units, A smooth surface is characterized by specular or 
mirror-like reflection with the angle of incidence determining the orientation of 
the reradiation pattern (Figure lA). If surface irregularities are a significant 
fraction of a wavelength, more energy will be scattered at angles other than near 
specular (Figure IS), If radar wave penetration and volume scattering is involved, 
sub-surface geometry also governs the reradiation pattern, A change in both the 
shape and relative magnitude of the reradiation pattern can be induced by increas- 
ing the dielectric constant. 

lurdien [1966, 1971] examined in the laboratory the effect of soil moisture 
on radar return, while MacDonald and Waite [1971] using commercially available 
imaging radars (K“bond multipolarization) showed that differences in gross soil 
moisture content could be determined qualitatively. Measuring the microwave 
return under natural conditions and with quantitative soil moisture and configuration 
information was begun in August, 1972 at the University of Kansas Center for Research 
Inc, Preliminary results of backscatter measurements [Ulaby, 1974] taken on two 
fields by a truck-mounted 4-8 GHz, FM-CW radar system (Table 1) indicated that 
radar response to soil moisture content is highly dependent upon the surface rough- 
ness, microwave frequency, and look angle. Ulaby et al, [1974] in a discussion 
concerning a new model based on the skin depth (attenuation) concept for soil 
moisture using the 4-8 GHz system, concluded thot within 10*^ of nadir, a strong 
and essentially linear relation was observed at all frequencies between the scatter- 
ing coefficient and the effective attenuation coefficient. They also noted that at 






Type; FM“GA^ 

ModuiaHng Wave Form; Triangular 

Center Frequencies; 4.7 GHz, 5,9 GHz, 7.1 GHz 

Bandwidth; AF 1 .2 GHz 

Transmitter Power; 5 watts 

IF Frequency; Fjp 87 KHz 

IF Bandwidth; AFjp 5 KHz 

Antennas; 

Height above ground 67 feat 

Transmitting antenna diameter 2.5 feet 

Receiving antenna diameter 3,0 feet 

Feeds ridged waveguide, dual polarized 

Beamwidths of the patterns product 

(Gj (©, $) . (0, 'P)) Elevation; {over 4-8 GHz) 

Azimuth; 3, 8^-2, 9° (over 4-8 GHz) 
0° “ 80° from the vertical. 


Look Angle Range; 0 




larger angles the sensitivity to changes in the effective attenuation coefficient and 
therefore in moisture content, was greater at the lower frequency. The skin depth 
ot the highest moisture values for the 7o1 GHz frequency is approximately one 
centimeter [Ulaby, etaU, 1974] which meant; 

1 « These measurements are less valuable for determining moisture values to 
significant depths. 

2, The validity of the "ground "truth” data is questionable due to the thin 
moisture region and the difficulty Involved in defining a rough surface. 

They concluded that lower frequencies and therefore greater skin depths appear 
necessary if the total moisture content for a given soil thickness is to be directly 
measured rather than inferred from the moisture in very thin surficial layers [Ulaby, 
et al., 1974]. Table 2 illustrates skin depth as a function of soil moisture and 
frequency. 

The measurement of the complex dielectric constant affords an excellent 
estimate of the soil moisture content due to its dependency to the content of free 
water in the soil. Moreover, at the lower microwove frequencies, Lundien [1966] 
has shown that the effects of soil type on the value of the dielectric constants are 
overriden by the effects of the free water content in the soil . Cihlar and Ulaby 
[1974] have noted that the sole parameter affecting the frequency dependence of 
the dielectric properties of moist soil is water. For a smooth surface, the power 
reflection coefficient is directly related to this parameter. However, as the surface 
roughness increases this optimum relationship is no longer true and a correction 
factor must be applied. Experiments by Batlivala and Ulaby [1975] have resulted 
in a detailed analysis of the effect of roughness on cP for bare fields. Their 
objective was to determine an optimum set of radar parameters (frequency, poiari" 
zation, incidence angle) such that the a° value of the ground is virtually indepen- 
dent of surface roughness but not at the expense of a reduction in sensitivity to soil 
moisture variations. The following discussion is a summary of their analysis. 

Experimentation was carried out using a truck-mounted 2 GHz - 8 GHz, 
microwave active spectrometer (MAS) system [Oberg and Ulaby, 1974] on three 
fields of considerably different surface roughnesses, but each denuded of any 
vegetative cover. The roughness was calculated using rms height with respect to 
the mean surface of the ground. The smooth, moderate, and rough fields were found 
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TABLE 2 Skin Depth 6 in cm for the Soil Moisture Profiles of Figure 5. 


Incidence Angle 


Number^ 

Date of 
Measurement 

Moisture 
Content ^ 

Frequency 
in GHz 

0° 

10° 

20° 

30° 

40° 

50° 

60° 

70° 




4.7 

6.2 

6.2 

6.2 

6.2 

6.2 

6.2 

6.2 

6.2 

1 

8/18/72 

4.8 

5.9 

5.8 

5.8 

5.8 

5.8 

5.8 

5.7 

5.7 

5.7 




7.1 

5.5 

5.5 

5.4 

5.4 

5,4 

5.4 

5.4 

5.4 




4.7 ■ 

3.3 

3.3 

3.3 

3.3 

3.3 

3.3 

3.3 

3.3 

2 

8/15/72 

15.8 

5.9 

3.1 

3.1 

3.1 

3.1 

3.0 

3.0 

3.0 

3.0 




7.1 

2.8 

2.8 

2.8 

2.8 

2.8 

2.8 

2.8 

2.8 




4.7 

1.4 

1.4 

1.4 

1.4 

1.4 

1.4 

1.4 

1.4 

3 

8/29/72 

24.0 

5,9 

1.2 

1.2 

1.1 

1.1 

1,1 

1.1 

1.1 

1.1 




7.1 

1.0 

1.0 

1.0 

1.0 

1.0 

0.9 

0.9 

0.9 




4.7 

U4 

1.4 

1.4 

1.4 

1.4 

1.4 

1.4 

1.4 

4 

9/5/72 

30,2 

5.9 

0.9 

0.9 

0.9 

0.9 

0.9 

0.8 

0.8 

0.8 




7.1 

0.7 

0.7 

0.7 

0.7 

0.7 

0.7 

0.7 

0.7 




4.7 

1.1 

1.1 

1.1 

1.1 

1.1 

1.1 

1.0 

1.0 

5 

8/25/72 

20.0 

5.9 

0,9 

0.9 

0.9 

0,9 

0,8 

0,8 

0.8 

0.8 




7,1 

0.7 

0.7 

0.7 

0.7 

0.7 

0.7 

0.7 

0.7 




4.7 

1.5 

1.5 

1.5 

1.5 

1.5 

1.5 

1.4 

1.4 

6 

c 

20.0 

5.9 

1.2 

1.2 

1.2 

1.2 

1.2 

1.2 

1.2 

1.1 




7.1 

1.0 

1.0 

1.0 

1.0 

1.0 

1,0 

1.0 

0.9 


^^Average moisture content of the top 5 cm of soil m^^ in percent by weights 

L 

Moisture profiles are shown in Figure 5. (After Ulaby, et al., 1974). 

^^Hypothetical case, moisture profile assumed constant with depth. 



to have values of 0,88 cm, 2,6 cm and 4.3 cm respectively. The three fields were 
chosen such that they most closely represented typical plowed fields for the rough 
category, and plowed“then"disced for the moderately rough category. The smooth 
field is a somewhat rarely occuring category but elimination from consideration 
would have implied a bias toward rougher fields. 

Each radar data set utilized 5 angles of incidence from 0° (nadir) through 
40° in 10° intervals for both HH and W polarizations at eight central frequencies 
between 2 GHz and 8 GHz (Table 3), It was necessary to use spacial and frequency 
averaging to reduce signal fading. 

Soil samples were obtained at 8 different locations for each field per data 
set at five depths; 0”1 cm, T“2cm, 2"5 cm, 5“9 cm, and 9“15 cm. Moisture con- 
tents were determined by the weight and by the bulk density methods. 

Surface roughness, as previously discussed, is a geometric property of the 
terrain and is expressed relative to the wavelength. For a field of a given surface 
configuration, an increase in frequency would imply a rougher field in terms of the 
electromagnetic rodiation impinging upon it. Figure 2 illustrates this point. It 
should be further noted that where the curves intersect, the effect of surface rough- 
ness is minimal; those values are 4° for 2,75 GHz, 10° for 5,25 GHz and 20° for 
7.25 <5Hz, 

When considering the effect of roughness on the spectral and the moisture 
responses of 0 °, the following conclusions are apparent; at nadir (Figure 3A) ct° 
displays little variation between 5,25 and 7,25 GHz but exhibits a decreasing 
response with frequency between 2,25 and 5,25 GHz for moderate and rough 
fields. At 10°, 0 ° for the rough field decreases with frequency but at a slower 
rate than for the same bandwidth at nadir, the moderately rough field displays an 
independence to frequency and the smooth field reverses slope,, while at 20° 

(Figure 3C), cr° for the smooth field increases rapidly. The tr° response to soil 
moisture is illustrated in Figure 4 (a,b)for three surface roughness profiles at 2,75 
GHz for two angles of incidence; 0° (nadir) and 10°. Slope is an indicator of the 
0 "° response to moisture and is referred to as sensitivity given in dB/.Ol g/cm , 

For example, the sensitivity of the smooth field (Figure 4a) is the greatest at nadir 
while the sensitivity of the rough field is the least. 

At 10° (Figure 4b) just the opposite observation holds. This would imply 
that the roughness effect can be minimized by operating at an angle between 0° 
and 10°, 
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TABLE 3. MAS 2-8 Sysfem Specifications 


Types 

Modulating Waveform; 

Center Frequencies: 

FM Sweep: AF 
Transmitter Power; 

IF Frequency; Fjp 
IF Bandwidth: AFjp 
Antennas; 

Height above ground: 
Transmitting antenna diameter; 
Receiving antenna diameter: 
Feeds; 

Effective Two"way beamwidth; 
Incidence angle range: 

Polarization: 

Calibration; 

Internal 

External 


FM-CW 

Triangular 

2.75, 3.25, 4.75, 5.25, 5.75, 6.25, 

6.75, 7.25 GHz 
450 MHz 

40 mW 
50 KHz 
6 KHz 

20 m 
91 ,5 cm 
91 ,5 cm 
Log periodic 

5.4° at 2.75 GHz/ 2.2° at 7.25 GHz 
0° (nadir)-80° 

Horizontal transmit“Horizontal receive (HH) 
Vertical transmifVertica! receive (W) 

Delay line 
Luneberg lens 
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Frequency (GHzl; 2.15 
Polarization: VV 


Frequency (GH’l; S. 25 
Polarization- W 



_i j 1 

20 3Q 40 

Angle of Incidence (Degrees) 

la) 2. 75 GHz 


Frequency (GHz): 7,25 
Polarization: W 



Angle of Incidence (Degrees) 

(bl 5. 25 GHz. 

Figure 2. Angular response of the scat" 
taring coefficient for the smooth, medium 
rough, and rough fields for high levels 
of moisture content at (a) 2.75 GHz, 

(h) 5.25 GHz, and (c) 7.25 GHz. (From 
Batlivala and Ulaby, 1975), 
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10 20 3 0 40 
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Scnllcring Coeflicicnt IdOl ScaUerlng C«;ff[cient (dB) 


Poiarizaiiott: vv 

Polarization: VV 

Angie of Incidence; 0® 

Angle of Incidence: 10® 


fiMS Height 
in cm 


RMS Height Soil Moi stu re 
In cm Ing/cm^lni;) 


?.25 Z.75 3.S 3.75 425 475 5,25 5.75 6.25 6.75 7.25 
Frequency tCHzJ 
(a) 0*’ (nadir! 


125 2L 75 3,25 3,75 425475 5.25 5.75 6,25 6,75 7.25 
. Frequency (GHzl 

it)} 10° Angle of Incidence 


Polarizations W 
Angle of Incidence* 2(f’ 



RMS Height 

Soif Moisture 


in cm 

iny/cm^lm;l 


0.SS 

0.333 


2.6 

0 350 


43 

0.360 





2^25 2.75 3.25 3.75 425 4 75 5.25 5.75 6.25 6.75 
frequency iGHz) 

(c) 20® Annie of Incidence 


Figure 3, Spectral response of the scatter" 
ing coefficient for smooth, medium rough, 
and rough fields for high level of moisture 
conteti|’ at (a) 0 (nadir), (b) 1G° and 
(c) 20 angle of incidence. (From Botlivala 
and Ulaby, 1975). 













Frequency (GHzl; 2.75 
Angle of Incidence: 0° 

PolarizalioH: VV 

• Smooth Surface (RMS Height * 0.88 cml 
▼ Medium Rough Surface (RMS Height = 2.6 cm) 
« Rough Surface (RMS Height = 4.3 cm) 

S * Soii Moisture Sensitivity (dB/0.01 g/cm^l 
R = Correlation Coefficient 



Frequency (GHz): 2.75 
Angle of Incidence; 10“ 

Polarization: VV 

• Smooth Surface (RMS Height = 0.88 cm) 

▼ Medium Rough Surface (RMS Height ■ 2, 6 cm) 
■ Rough Surface (RMS Height = 4.3 cm) 

S ■ Soii Moisture Sensitivity (dB/0. 01 g^cm^) 

R = Correlation Coefficient 



(b) 10“ Angie of incidence 


Figure 4. Soil moisture responses for the three surface roughness profiles at 2,75 GHz for (a) 0° (nadir) 
□nd (b) 10 angle of incidence. (From Batlivala and Ulaby, 1975). 



The effect of surface roughness on fhe choice of optimum system parameters 
is seen in Figures 5 " 7 which depict calculated linear regression lines of as a 
function of soil moisture content. As the frequericy is increased theiradar response 
is more sensitive to roughness at a! I soil moisture conditions but a decrease in sen*“ 
sitivity occurs for the soil moisture, ThuS/ operating at 0° (nadir) is undesirable 
due to this lai^e change in <P with respect to rms height especiolly at 0,35 g/cm^. 
(Note: 14 dB at 2.75 GHz, 18 dB at 4.75 GHz, and 16 dB at 7.25 GHz). 

At an incidence angle of 10°, the response to the variation in surface 
roughness (Figure 6) between cr° for the smooth field and cP for the rough field 
decreases as the frequency is increased from 2.75 GHz to 7,25 GHz. The optimum 
condition is one where the moisture curves for cP are horizontal lines, parallel, 
and at maximum separation. The frequency which best satisfies these three condi- 
tions is 4,75 GHz, 

o o 

Variations in cr at 20 (Figure 7) also decrease as frequency increases, but 
the data indicates that for roughness to become an independent variable, a frequency 
greater than 7,25 GHz is necessary. Increasing frequency or angle results in a 
reduction in the magnitude of sensitivity, the latter represented in terms of separo- 
tion between the curves. Thus it appears that the optimum frequency is 4.75 GHz, 
W polarization at 10° Incidence angle. 

To expand the number of available angles between 0° and 20°, a non-linear, 
interpolation technique was employed to evaluate the cr° values at 3.3°, 6,7°, 

13.3° and 16.7°. This procedure was applied for each data set at each frequency 
and polarization. Hence a new data base consists of 8 frequencies, 9 angles and 
two polarizations? a linear regression analysis was performed for each data set. 
Optimum correlation coefficient, sensitivity and frequency plotted as a function of 
incidence angle for moderately rough and rough fields are graphically illustrated 
in Figure 8, The optimum correlation is at the maximum (0,89) at 10° and at 2,75 
GHz, for both polarizations. The optimum frequency in the 0° to 20° angular 
range varies between 2.75 GHz and 3,25 GHz, When the data is combined for 
the moisture response for smooth, moderate, and rough fields (Figure 9), the 
optimum parameters exhibit a strong dependence on incidence angle. 

These conclusions were based upon optimum sensitivity and correlation coef- 
ficient considerations? Figures 8 and 9 indicate an angulor range of approximately 
7° to 15° as optimum. An optimum frequency of about 3 GHz is suggested by 
Figure 8 while Figure 9 suggests a frequency of about 5 GHz; 4 GHz is recommend- 
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Angle of Incidence: 0® 
Freguenty tCHtf: 2.75 
Polarizalion: VV. 

Soil Moisture 



• 0.25 

0.35 



_) t ! ! 

2.0 3.0 40 5.0 

RMS Height (cml 

l3) 2. 73 GHz, nadir 


Angie ol Incidence: 0® 
Frequency iGHzh ^.75 
Polsrizalion: VV 
Soil Moisture (ms 1 gfcm®! 
■ ao5 



Angle of incidence: 0® 
Frequency IGHzh 7.25 
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Figure 5, Scattering coefficient as a 
function of surface roughness at nadir 
for four moisture conditions at (a) 2*75 
GHz, (b) 4.75 GHz, and (c) 7.25 GHz 
(From Batlivala and Ulaby, 1975). 
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Figure 6, Scafi'ering coefficient as a function of surface roughness at an angle of 
ihciiience of 10° for four moisture conditions at (a) 2,75 GHz, (b) 3,25 
GHz, (c) 4,75 GHz, and (d) 7.25 GHz, (From Batlivala and Ulaby, 1975). 
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Figure 7, Scattering coefficient as a 
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GHz, and (c) 7,25 GHz. (From Batlivala 
and Ulaby, 1975). 
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ed as a compromise. Furthermore, both HH and W polarizations displayed 
approximately the same moisture response. 

The preceding discussion was covered in considerable detail because it 
represents the most recent work in the field of soil moisture detection utilizing 
radar in which ground data was an integral part of the analysis. Many of the 
following recommendations for monitoring soil moisture from space are the 
direct result of this investigation. 

One further point should be noted: while the above results are applicable 
for fields void of any vegetative cover, similar observations have been obtoined 
when considering vegetation [Ulaby, 1975; Ulaby et al,, 19751. The correlation 
coefficient between cr° and soil moisture content is greatest at a 10*^ angle of 
incidence at 4.7 GHz even though the presence of vegetated cover tends to 
reduce the sensitivity of cP to soil moisture variations, 

3,0 CONCLUSIONS 


A tool of great potential exists whereby soil moisture conditions could be *> 
monitored on a regional scale if the proper airborne platform were utilized. If 
the objective were to determine a temporal change in soil moisture through repeti“ 
tive coverage of a given area, a revisit time of 1 to 30 days, the latter being a 
maximum, is required of any water resources satellite program (depending upon the 
area, parameters interested in, etc.) (periodicity of coverage obtained from tables 
of measurable parameters from a report entitled, "Satellite Data Collection User 
Requirements: Report of the NASA Ad Hoc Committee on Hydrology). A short 
revisit time is necessary due to the extreme changability of this surface phenomenon. 
To realize the maximum useage of data retrieved in improving hydrological fore- 
casting and water resource management, an average revisit time of 4 days is sug- 
gested, This requirement leads to the necessity of using very large swath widths, 
about 400 km, to ensure full coverage of the continental United States plus 
Alaska, Based upon the optimum angular range of incidence angles (7° - 15°) 
for monitoring soil moisture conditions, an orbiting spacecraft at 1 100 km would 
produce a ground swath width of about 160 km. It should be noted that design 
considerations dictate an image swath width based upon an angle of incidence 
range of 7° to 22°. While the overall dimensions of the swath approximates 300 km, 
any point beyond the optimum angle (beyond 15°) would fall into the ambiguity 
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category. Exactly where the cutoff is for useable data beyond the 15*^ angular 

i I 

range is speculative. This problem can only be resolved when more data becomes 
available. 

The above discussion centers on the criteria of repetitive coverage in order 
to insure accurate monitoring of soil moisture variations on a temporal basis. If ] 

on the other hand, the objective was to map soil moisture content from a sirgle 
mission, the complications induced by the 4 day revisit requirement are eliminated. 

However, the optimum angular range of incidence angles necessary for maximum I 

data retrieval remains the same. Moreover, image swath width is still a major 
consideration from the standpoint of mission planning since.it may now be necessary 
to fly more than one pass for a given area to obtain complete coverage. Such a 
consideration would be unfeasible for same*“day coverage if an orbital platform were 
utilized. 

The proposed system resolution of 100 meters [Claassen, 1975] should pro- 
vide sufficient information for both regional analyses and some local detail. 

Greater spacial resolution is impractical due to power requirements [McMillan, ; 

1975]; poorer resolution would hamper data analyses for all but gross variations 
in soil moisture. 

In conclusion, the following recommendations are submitted for application 
toward an operational system. The recommended radar parometers are; 

7^ - 15° 

4 GHz 
HH or W 
100 meters or better 

A revisit time of 4 days is recommended. 


optimum angle of incidence range; 
optimum frequency; 
polarization; 
resolution; 
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ABSTRACT 


Improved snowpack melt forecasts are necessary to ensure the success- 
ful management of any water resource program where the major source of 
water is derived from snow melt runoff. Current technology has provided 
techniques whereby snowpack monitoring is feasible from spaceborne, remote 
sensing systems. ERTS has yielded improved methods for snowpack, data 
acquisition. Similar results have been obtained utilizing the NOAA-VHRR 
and the ESSA-9 systems. However, problems of resolution, cloud cover 
restrict ions, canopy problems, and non-con tinuous snow cover ambiguities 
have often rendered the images less useful than originally expected. 

Experiments with active microwave sensing systems have been performed; 
the University of Kansas 1-8 GHz Microwave Active Spectrometer (MAS 1-8) 
system was used to measure the backscatter response of snow covered ground 
between 21 February and 23 April 1975. Simultaneous ground truth data con- 
sisted of soil moisture, soil temperature profile, snow depth, snow temp- 
erature profile, and snow water equivalent. The preliminary results of 
the experiment indicate that the key parameter effecting the radar return 
is snow wetness when density stratification is not present within the 
snowpack under investigation. 
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RADAR MEASUREMENT OF SNOW — 

STATE OF THE ART 
Bradford C. Hanson 

Snowpack water resources are of considerable importance in the man- 
agement of irrigation systems, flood warning systems, hydroelectric- 
power schedules, and municipal and regional water supplies. Snowpack 
melt often provides the major source of water for power and human consump- 
tion for many areas throughout the world. In California for example, the 
Sierra Nevada snowpack accounts for more than 1/2 the total water supply 
and about 1/3 of the electric energy. Monitoring the seasonal snowpack 
parameters of surface temperature, total thickness, water equivalent, 
liquid water content, plus their continuous changes with respect to each 
other is necessary to ensure the successful management of any snowpack, 
water resource program. 

An accurate water supply forecast depends upon an adequate and timely 
input of all necessary variables. The present technique employed for 
snowpack data acquisition is snow surveys. This traditional method con- 
sists of establishing a number of sample points within a given watershed 
for which physical sampling is performed. These measurements are then 
related to runoff by one of several methods, the historical normal, an 
index method, a water balance method, or by hydrologic methods. This 
information from aggregate spot locations has been and continues to be 
the basic snowpack evaluation technique for synthesizing the snow cover 
and resultant runoff for an entire watershed. However, too few data 
points are normally available to make reliable estimates of the snow 
cover and its physical condition. Adding additional ground stations Is 
a costly venture. In mountainous regions of the U.S. for example, mil- 
lions of dollars are spent each year at fixed locations to measure the 

snowpack. Moreover, improved forecasts are estimated to be worth 10' to 
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10 dollars per year to water users in the western United States alone 
(Useful Application of Earth-Oriented Satellites, 1969). Another problem 
is the timely input of data at regular intervals plus the necessity of 
data upon demand. Watershed areas in the mountainous regions of Arizona 
often undergo rapid change within a time span of a few days to a few hours. 


Thus one must not only monitor vast areas incorporating an infinite number 
of sample points, but must also monitor on a periodic basis. The problems 
encountered from such a technique are numerous, but until recently, other 
methods were impossible due to the lack of the necessary technological 
advances. 

Current technology has introduced the feasibility of airborne and 
spaceborne remote sensing systems for monitoring the snowpack. Water- 
shed snowpack and streamflow data obtained and transmitted by ERTS have 
been used by Warskow et al., (1975) in operational and water management 
decisions for the Salt River project. However the satellite image gave 
only \k% reliability and was therefore used only to provide an overview 
for establishing snow course sites. Other problems have been encountered 
when applying ERTS data to snowpack prediction techniques; cloud cover is 
usually the rule rather than the exception especially immediately pre- 
ceding and following a major snow storm. Bare rock often appears as 
snow cover (Washichek and Mikesell, 1975), canopy problems (Limpert, 1975; 
Kathibah, 1975), and Imagery density and contrast Variations which are 
not due to ground data (Aul and Ffolliott, 1975) are also encountered. 
Moreover at optical wavelengths, information retrieved is only from the 
surface of the snowpack; its physical condition at depth can only be 
inferred. 

The National Environmental Satellite Service (NESS) began producing 
river basin snow maps on an experimental basis in the spring of 1973* The 
primary source of data was from the VHRR (Very High Resolution Radiometer) 
aboard the NOAA polar-orbi ting satellites at 0.55 to 0.75 pm to 10.5 to 
12,5 urn. Fog, cloud cover, and severe resolution restrictions rendered 
the data less useful than originally expected (Schneider, 1975). The 
U.S. weather service in Alaska employed NOAA-VHRR enhanced IR satellite 
imagery as an aid in the acquisition of daily synoptic snowmelt informa- 
tion for use with river forecasts (Seifert, et al., 1975). Again the 
major problems were cloud cover and image resolution, the latter being 
one mile which made variations within the snow cover impossible to detect. 
Resolution for the ESSA-9 system has been reported to be h kilometers 
(Rango, personal communication). Other difficulties encountered include 
the ambiguity between snow and clouds, identifying snow In areas of 
mountain shadows, and mapping snow where the snowcover is not continuous 
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(Barnes and Sitiallwood, 1975). 

Edgerton et al., (1971) recorded microwave measurements of snow at 
1.^3s ^.39, 13.^ and 37 GHz and reported a direct relationship between 
brightness temperature and the water equivalent of dry snow (Figure 1). 

A general decrease in brightness temperature as total snowmass increased 
was also noted. This trend was confirmed in a later experiment in which 
the water equivalent was found to be the main factor influencing the 
brightness temperature from snow. 

Cosgriff et al . , (i960) completed one of the original measurements 
of snow cover with a truck mounted CW-Doppler scatterometer operating at 
X-band and K-band frequencies. Their results indicated that in general, 
spow had the effect of covering the backscatter from the terrain (Figure 
2), Later Moore (1972) plotted this same data versus water content and 
illustrated a 1 inear relationship between the normalized scattering coef- 
ficient and total water content (Figure 3). These data however, were at 
various temperatures and had insufficient ground truth to determine whether 
the results were actually due to the free water content within the snow- 
pack or due to the underlying soil moisture. Vickers and Rose (1972) 
utilized a 2.7 GHz short pulse radar to investigate snowpack stratification. 
Their method has the drawback of needing ground based calibration for 
either the snow density or snow depth for the other variable to be remote- 
ly determined. Linlor (197^) -measured the microwave reflection from various 
snow depths and found considerable variation with depth (Figure 4) but made 
no attempt to analyze the response. Waite and MacDonald (1970) observed 
differences in signal return between fresh snow and older perennial snow. 
They noted that in the past measured differences had been primarily attri- 
buted to moisture content variations within the snow pack, but in view of 
their data, they felt that these differences were probably due to volume 
scattering from inhomogenei ties within the perennial snow. Moreover, 
cross polarization (HV) tended to bring out this effect more than like 

V 

Ik 

polarization (HH) . Unfortunately, ground data did not exist for further 
confirmation of their results. The lack of simultaneous ground data in 
conjunction with aerial measurements over the snow pack Is the major prob- 
lem with most of the snow cover experiments to date. One notable exception 
is data acquired during the winter of 1974 and 1975 at The University of 
Kansas. 
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Figure 1. Brightness ’’emperature Versus Water Equivalent Dry Snow, 

Crater Lake, Oregon, 22 March 1970 (Edgerton, et al . , 1971). 
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Figure 2. Effects of Snow on Gamma at Ku-Band 
(Cosgriff, et al., I960). 
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Radar backscatter data were acquired with a truck^mounted 1-8 GHz 
Microwave Active Spectrometer (HAS 1-8) at eight frequencies within the 
1-8 GHz range, for three polarizations (HH, HV, VV), and for incidence 
angles between 0® and 70°. Table 1 summarizes the system specifications. 

Ground data (Hanson, 1975) was acquired in conjunction with the 
spectrometer measurements. Data collected consisted of soil moisture, 
soil temperature, air temperature, snow depth, snow temperature, and the 
stratigraphic profile within the snow pack. Since the experiment was 
conducted during the winter dormancy period, the height and physical 
condition of the ground cover was essentially constant. 

Figures 5 and 6 illustrate scattering cross section, spectral responses 
for three distinct ground conditions. The maximum accumulation of dry snow 
encountered during the experiment was 15 cm; the effect of this minimal 
snow depth is demonstrated by Figure 5 which compares the response for 
Run 4 (no snow) and Run 6 (snow cover - 15 cm dry powder) . This is con- 
sistant with preliminary results obtained by Stiles and others (1976) 

(Table 2) which indicate that only a small effect should be observed for 
dry snow at depths of 15 cm or less since the skin depth for snow is great- 
er than the overall snow depth. 

Immediately before Run 8, a fine sleet Intermixed with a light mist 
accumulated on the surface of the snow. Although snow wetness could not 
be measured, the surface was damp. Figure 6 illustrates the cr° response 
to the control set, Run 4, and the wet snow cover set. Run 8. The skin 
depth of "wet" snow at 1.0 GHz is 78 cm (Linlor, 1972; Table 2). Thus at 
1.2 GHz, the 12 cm of wet snow (Run 8) in Figure 6 acts primarily as an 
attenuator of the backscatter from the underlying wet soil (Stiles et al., 
1376). The skin depth for the wet snow cover at 7.25 GHz is comparable 
to or smaller than the 12 cm depth (Table 2) and therefore it is doubtful 
that there is any significant contribution to the 7-25 GHz response of 
Run 8 due to the underlying soil (Stiles et al., 1976). 

The major conclusion derived from this experiment is that the key 
parameter effecting the radar return Is snow wetness (Stiles et al., 1976). 
However it should be noted that snow depths were minimal and that cMmatic 
variations were such that density changes within the snow cover were not 
observed. In areas where greater snow depths occur for longer periods of 
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Table 1. HAS 1-8 Specifications 


Type: 

Modulating Waveform: 
Center frequencies: 


FM sweep: AF 

Transmit power: 

IF bandwidth: AF,_ 

I F 

Antennas: 

Height above ground: 
Transmit antenna diameter: 
Receive antenna diameter: 
Feeds : 

Effective system beamwidth: 
Incidence angle range: 
Polarization: 


Cal ibration: 
Internal 
External 


FM-CW 

Triangular 

I. 2, 1.6, 2.25, 3.25, ^.25, 5.25, 
6.25, 7.25 GHz 

400 MHz 
loo mW 
49.5 KHz 

II. 5 KHz 


19 m 
1 .22 m 
1 .22 m 

Dual polarized switchable 
1-12 GHz log periodics 
12“ at 1.2 GHz to 2“ at 7.25 GHz 
0“ (nadir) - 80“ 

Horizontal transmlt-horizontal 
receive (HH) 

Horizontal transmit-vertical 
receive (HV) 

Vertical transmit-vertical 
receive (VV) 

Coaxial delay line 
Luneberg lens 
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Table 2. Skin Depths for Snow 


Description Frequency (GHz) 

iL 


tan 6 

Skin Depth 

Source 

Dry Snow 

37 

1.9 

.05 

— 

7 cm 

Edgerton 

Dry Snow 

13.6 

2.76 

.03 

— 

38.8 cm 

Edgerton 

Wet Snow 

1 .0 

1.5 

1 

10 -' 

78 cm 

Linlor 

Wet Snow 

8.0 

1.5 

— 

10 -' 

9.7 cm 

Linlor 

Dry Snow 

1.0 

1.32 

.0002 

— 

548 m 

Evans 

Dry Snow 

8.0 

1.32 

,0002 

— 

68 m 

Evans 

Foam 8.5^ H^O 

1 .83 

— 

— 

— 

43 cm 

Linlor 

Foam 8.5% H^O 

8.0 

— 

— 

— 

3 cm 

Linlor 

Saturated Snow 

10 

— 

— 

— 

.9 cm 

Linlor 

Freshly Fallen 
Snow 

3 

1.2 

— 

.00029 

100 m 

Vickers 

Freshly Fallen 
Snow 

10 

1.26 

— 

. 00042 

20.2 m 

Vickers 

Hard Packed Snow 

3 

1.50 

— 

.0009 

28.9 m 

Vickers 

Ice 

3 

3.2 


.0009 

19.8 m 

Vickers 


(After Stiles et al., 1976). 
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time and stratification is evident, the key parameter effecting radar 
return may no longer be snow wetness. 

In conclusion, the technology exists whereby orbital platforms 
containing accurate sensor packages are feasible for measuring snow pack 
conditions. The major stumbling block at the present time is the lack of 
specific information regarding the proper system design such that optimum 
data retrieval can be realized. This much-needed information can only 
be obtained after further experiments. 
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ABSTRACT 


Flood plain and flood inundation mapping, surface water mapping including 
estuary and wetland surveys, and flood area assessment has received considerable 
attention during the past decade because of man's increasing environmental aware" 
ness. Great strides have been achieved in creating new, more efficient methods for 
monitoring these parameters; most -^aparent are the innovations resulting from the ERTS 
and the SICYLAB programs. Active microwave sensors possess an even greater potential 
due to their independence of the operational constraints inherent to visible light“IR 
systems and their resolution superiority over passive microwave systems. An orbital 
platform containing a synthetic-aperture radar system would provide ah additional, 
extremely valuable, tool for any water resource mapping program. 

Optimum system parameters required for maximum data retrieval are difficult 
to establish due to an inadequate research data base. Based on radar theory and 
existing limited research, the following recommendations ore submitted as realistic 
sensor instrumentation requirements: 


frequency: 

X43and 

polarization: 

HH 

depression angle 


limitation: 

less than 75* 

resolution: 

12 meters 

revisit time: 

6 days 



1.0 INTRODUCTION 


Water percolates downward, displacing interstitiai air within the zone of 
aeration until ail pores are permeated. With additional moisture, downward infiltra- 
tion will be impeded if the soil is in a state of maximum saturation. Additional surface 
water under these conditions exists either at the alr-soM interface as free-standing 
water, or becomes runoff. Excessive amounts of the latter often cannot be contained 
within the confines of the river channel; when a drainage system can no longer 
accomodate the increased volume, flooding occurs. Floods are an extreme case of 
surface water runoff, but for any water resources management program, the ability 
to accurately Identify and map flooded areas, and delineate those areas of potential 
flooding is of prime importance. 

Most surface-water runoff is contained within a given drainage systen or is 
entrapped by an inland water body, the latter comprising reservoirs, lakes, ponds, 
etc. Surveys involving the measurement of the total area of surface water over wide 
regions, continental or global, would be of scientific as we!! os practical use. The 
amount of water subject to evaporation, populoti on -center accessibility to a water 
supply, more accurate indices of the total water available for consumption and ifs 

change with respect to time, would be more rapidly obtained utilizing remote sensing 
techniques. Volume information is of equal importance. Volumetric data can be 
utilized in the analysis of a local/regional water table, a parameter indicative of 
the water storage capacity for a particular area. Thus an equally important factor 
in a water resources management program is the ability to chart and measure surface 
water (lakes, ponds, marshes, etc.) on a local and on a regional basis and at 
regular intervals. The latter stipulation is necessary if monitoring changes with 
respect to time are to be realized. 

The purpose of this paper is twofold: (1) to explore the potential use of an 
active microwave sensor to delineate various surface water phenomena and (2) to 
identify the proper radar design parameters for a system operating from an orbital 
platform. 

2.0 PREVIOUS INVESTIGATIONS 


Many investigations [Anderson, et a!., 1973; Klemas/ etol,, 1973; Mairs, et ai. 
1973; Reeve:?, 1973] have dealt with flood plain and flood inundation mapping, surface 



water mapping and flood area assessment including estuary and wetland surveys utiliz" 
ing ERTS"1 data. Extensive work has been performed using ERTS because the Images 
are readily attainable and because open areas of surface water are easily delineated 
by their strong contrast with the surrounding land mass at 0,8 - 1 ,1 ym. The major 
drawbacks have been the rather long repeat cycles of 18 days (? days now possible 
with two satellites), resolution restrictions [Reeves, 1973] and restrictions imposed 
by inclement weather conditions during daylight hours. Other problems arise when 
considerations of wetland mapping are involved. It is not always possible to delineate 
with precision the marsh-water interface and the upper wetland boundary though 
certain plant communities may be useful for this [Anderson, et al., 1973; Klemas, et 
al„ 1973;Mairs, etal., 1973], 

Radar imaging systems are not restricted by the problems previously discussed. 
Daylight operation is no longer a restriction with radar, nor does cloud cover impede 
data retrieval. Moreover, fine resolution is achievable from any altitude since 
resolution is independent of range for a synthetic-aperture radar. An added benefit 
is its ability to delineate physiographic characteristics which might otherwise go 
unnoticed on aerial photography (Figure 1). Note the channel configuration and the 
ossociated meander scars; surface features such as soil variations, agricultural patterns, 
vegetation, and land use boundaries are also vividly expressed. Note also that intra- 
channel physiographic features such as mudflats are discernable to a certain degree 
from the surrounding mass of water. 

Radar has been suggested as a useful tool for the detection and mapping of 
inland water bodies as far back as 1967 [McCoy, 1967] due to its characteristic 
appearance on imagery, Roswell [1969] concluded that lakes larger than 8 acres 
in area could be detected in well drained, lowland areas on AN/APQ-97 imagery; 
Simpson [1969] obtained similar results for an area in New England, 

Radar (X“band) has been utilized for flood warning and damage assessment 
on a very limited basis, Rydstrom [19701 presented several techniques for identifying 
flooded areas. He noted the strong contrasts in returns between flooded and non- 
flooded fields and that breaks in levees were identifiable due to a disruption in the 
high return generated from the levee. Dams and associated spillways displayed high . 
returns during normal pool, but during periods of high water when the spillway is 
active, no return is recorded for the spillway structure [Rydstrom, 1970], 
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FIGURE 1. AN/APQ-97 radar imagery — Mississippi River Valley near Vicksburg 
Mississippi . 









3,0 RADAR SIGNAL RETURN 


Depending on the relative roughness of the terrain, the radar return signal 
from vegetation, soil, rock, etc. is generally greater than the return signal from 
water. Radar signals are normally returned from the terrain to the receiver by a 
scatterh g e^ocess (reradiation) with the intensity of radar return (signal strength) from 
the temin Untermining the relative degree of brightness, on the radar image. Parameters 
that influence radar return are frequency, polarization, complex dielectric constant, 
incidence angle, and surface roughness. 

The local angle of incidence, 0, is the angle formed between an impinging 
beam of radar energy and a perpendicular to the imaged surface at the point of 
incidence (Figure 2A). The angle between a line from the transmitter to a point on 
the terrain, and a horizontal line passing through the transmitter is the depression 
angle (ct). On flat terrain there is a continuous change in the angle of incidence 
from a maximum at for range to a minimum, near normal incidence, in the near 
range. This relationship is modified with the introduction of a slope to the terrain 
surface (Figure 2B) . 

The energy incident on a horizontal terrain surface is "specularly" and/or 
"diffusely" reflected in varying proportions depending upon the roughness of the 
terrain. Surface roughness is a geometric property of the terrain expressed relative 
to wavelength. Surfaces with micro“relief much less than a quarter-wavelength appear 
smooth (no return) whereas surfaces with micro-relief on the order of a wavelength or 
more, appear "rough" . A smooth surface is characterized by specular or mirror-iike 
reflection with the angle of incidence determining the orientation of the rerodiation 
pattern (Figure 3A). Under these conditions, the reflection obeys Snell's Law (angle 
of incidence equals angle of reflection) with virtually all the reflected energy being 
contained within a small angular region about the "specular" angle. When water 
bodies such as lakes and rivers are imaged, they act as near-specular reflectors and 
as such direct most of the transmitted energy away from the receiver (Figure 4B) 
except neor vertical incidence (0 = 0°) where strong bockscotter is recorded (Figure 
4A). Numerous investigators have noted a strong return in the extreme near range 
(maximum depression - minimum incidence angle) for ocean surfaces (Figure 5A); 
however, most energy is reflected away from the radar for incidence angles exceeding 
10 to 15 degrees [Hanson and Dellwtg, 1973] . 
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FIGURE 2. THE RELATIONSHIP BETWEEN ANGLE OF INCIDENCE, 
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FIGURE 4. APPROXIMATE RERADIATICN PATTERN FOR 
REFLECTION FROM A SMOOTH SURFACE. 
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FIGURE 5. AN/APQ“97 radar imagery — Golden Gate area 
San Francisco, California. 
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The scatter of radar signals from disturfaed water surfaces is usually less than 
that from land, but can be a confusing factor in distinguishing land“water boundaries, 
Backscatter from waves is governed largely by the height of the wind^drlven capillary 
waves, with modulation by the tilting of these tiny wavelets by the larger waves. 

When the wind is high, the backscatter approaches that from land in strength. With- 
in 15® or so of the vertical^the role of the capillary waves diminishes and the return 
(much stronger than at larger angles) is determined by the larger waves. The typical 
variation of radar return in an image is indicated in Figure 6A, B, C, A ; in Figure 
6DA no return is observed in the near zone, probably because of a very calm surface. 

The strength of returns from water is important from the standpoint of surface 
water mappability. For example, area and volumetric information can only be real- 
ized if the shoreline for a particular water body can be delineated. Under certain 
circumstances, shoreline recognition and mapping would be impeded if the area under 
investigation was situated in the extreme near range; accuracy in shoreline location 
would be minimal due to a high return from both the land and the surface of the water 
(Figure 7; A). 

The discussion thus far has been concerned with high microwave frequencies, 
on the order of 10 GHz to 35 GHz, X-band to Ka-band, Existing limited research, 
compiled at these frequencies, suggests that any microwave frequency between these 
two end-points will suffice for detecting and mapping free-standing water (lake, 
reservoir, pond, etc*) within a system's resolution capability. Decreasing the fre- 
quency usually causes a surface to appear smoother; smoothing the effects of surface 
roughness instigates a decrease in the amount of energy backscattered by a target. 

For water the effect would be a decrease in the signal return at incidence angles greater 
than zero degrees even though the physical characteristics of the capillary waves 
remoined constant. Thus It may be possible to override the effects of increased radar 
return due to the capillary waves and hence may be possible to reduce the problem 
discussed for Figure 7. Unfortunately, comparative imagery, L-band versus a higher 
frequency system such as X-band, was not available for our study to either confirm or 
discount this contention. 

One study was undertaken to determine the feasibility of utilizing multi- 
plexed, synthetic-aperture X- and L"band radar to obtain specific data over varying 
terrains [Drake et al,, 1974] , For the purpose of water resources management, their 
analysis indicated that X“band imagery permits identification of small nonlinear and 
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FIGURE 6. AN/APQ-97 radar imagery Illustrating radar return variations for 
water bodies in the extreme near range. 
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narrow linear open wafer areas; L"band imagery e:>diibited a more subdued response 
to these features which often rendered the imagery useless for identification purposes. 
Analysis of shorelines indicated that radar at X“band and shorter wavelengths was 
superior to L"band imagery in every respect; shoreline delineation is ambiguous vdien 
the peripheral vegetation is low and of even height when imaged with an L"band 
system^ but is easily located on X-band imagery [Drake et al., 1974]. Pads of water 
lilies are faintly indicated on both X- and L”band imagery as well as hyacinths and 
to some degree reeds; differentiation is only possible by utilizing both polarizations 
(like and cross) [Drake et al. , 1974] . 

The radar return generated for flood inundated areas at differing frequencies 
should depend upon the surface characteristics of the water, the nature of any floating 
materials, and ,'hether all ground cover is completely submerged, or whether some 
is protruding above the surface of the water. At low frequencies, for example, radar 
return for inundated areas with considerable floating surface material (logs, trees, 
etc.) or an area in which surface water does not totally conceal all ground cover 
(similar to a marshy condition), may differ considerably from that generated by a 
higher frequency system such as Ko“band. Moreover, it has often been noted that 
islands could not be positively identified from masses of floating vegetation at either 
X“band or L-band frequencies [Drake et al., 1974]. 

The boundary between swamp or marsh and open water presents a different 
problem; Roswell [1969] observed that the boundary is often diffuse. He noted that 
the gradation of gray tones between open water, water/vegetation intermixing, and 
non“water surfaces creates difficulties when establishing definite boundaries. The 
area exhibiting high return (Figure 8A; Areas X) peripheral to the river are wetlands 
dominantly vegetated by the common reed, Phragmities communis , interspersed by 
Wild Rice , Zizania aquatica . The; vegetation dominant in the areas of moderate 
return characterized by the speckled texture (Figure 8A; Areas Y) is that of tall, 
slender Blue joint grass. Both vegetation habitats generally exhibit definite, recog- 
nizable boundaries relative to each other and to the land-water interface , An 
exception to the latter occurs where water and sparse vegetation are intermixed 
(Figure 8A; ^ ) thus creating a zone of transition. Definite boundaries are easily 
established for terrestrial vegetation habitats where plant diversity is prevalent. For 
example. Areas Z (Figure 8A) are dominated primarily by Black Spruce, Piceg 
mariana, of uniform height (approximately 6 meters) and density. These areas are 
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FIGURE 8. AN/APQ-97 radar imagery Illustrating return variations for 
marsh and wetland regions. 
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easily deiineai'ed from the surrounding vegefaHon (primarily deciduous forest) by 
{heir lower radar return intensities and from nearby water bodies by their textural 
characteristics. The same observation holds when considering marsh areas 
dominated by Black Spruce in other locales (Figure 8B; Area. Z). Moreover Area X 
(Figure 8B) is primarily Biuejoint grass interspersed with the common reed and inter- 
mixed by minor amounts of wild rice. Thus it appears that definite boundaries can be 
established in most instances, the exception being water intermixed by sparse 
vegetation. The effect of incidence angle on such discriminations is unknown. 

Observations by Drake and others [1974] indicate that strong returns from 
marsh areas on X"band imagery are reflections from the top of the vegetation. At 
L-band frequencies, penetration apparently occurs; marsh reeds up to 5 feet above 
the surface of the water are penetrated, resulting in specular reflection. One further 
observation was noted: marsh areas were often confused with certain types of range™ 
land and agricultural land on both L“band imagery and X-band imagery but especially 
the latter. It must be noted however, that the observations by Drake [1974] were 
performed on test sites located in Brevard County, Florida; vegetation differs 
dramatically with respect to its physical size, shape, and other characteristics depending 
upon region. Hence, vegetation characteristic to a marshland in Florida may not be 
the same as that which is characteristic of marshes in other parts of the country. Thus, 
observations for marshlands in Florida may or may not hold for marshes in northern 
Minnesota . 

The low return observed for Figures 8A and 8B, Area Z, is ihe result of the 
vegetation characteristic to these environmental niches; the Black Spruce of uniform 
height and density has created a smooth surface relative to the wavelength involved, 
hence the low return. The high return noted for Figures 8A and 8B, Area X and the 
moderate radar return observed for Area Y Is due to differences In vegetation. For 
example, the Reed grass, PhragmIties communis , are erect plants which vary from 2 
metere to 4 meters in height, possess flat blades 1 cm to 5 cm wide and 15 cm to 40 cm 
long, grow in dense clusters, and are often confused with wild rice (Fassett, 1940). 

Even though the water level was 2.7 feet above normal pool (U .S. Corps of Engineers, 
Grand Rapids, Minnesota) when the radar image was flown (October 15, 1965), the 
strong radar return (Figure 8A; Area X) is the result of backscatter from the plant 
community; these stout, erect plants were partially submerged but retained their density 
and physical characteristics above the surface of the water. The more moderate 
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refurn infensIHes from fhe Blue jo inf grass tfioughf fo be Calamagrosfis canadensis 
(Figure SA; Area Y), is also the resulf of backscafter from fhe vegeiafion. Bluejoinf* 
are fall, siender grasses, seldom jxceeding 2 meters in height, possess flat blades 
less than 1 cenfimefer wide and 20 cm to 40 cm long and often form dense ground 
covers (Fassett, 1940). When partially submerged however, the Bluejoinf grass does 
not produce a radar return signal as strong as that for the Phragmities community 
because the canopy of Bluejoinf grass above the surface of the water is not as dense 
and thereby probably permitting a small portion of the impinging radar signal to be 
reflected from the water surface . 

Frequency is a parameter which also appears to control the degree of back- 
scatter observed for different topographic situations. For marsh and wetland areas 
in Florida, L-band appears superior to X-band because of its penetration capability; 
for open water areas and shoreline delineation, X-band is a better choice 
[Drake et al., 1974]. For similar areas in northern Minnesota, microwave frequencies 
lower than X-band may create boundary problems. In fhe case of Area Z (Figure 8A), 
lowering the microwave frequency may result in an obscured boundary since longer 
wavelengths tend to smooth surface features; hence differences in radar return are 
more subdued. For areas where diffuse boundaries are a problem, lowering the 
frequency often results in a more clearcut situation; differences between the non-water 
surface and fhe intermixing zone become less evident and when compared to the 
open water zone a definite boundory becomes apparent [Roswell, 1969]. Thus there 
appears to be a tradeoff in retrievable data when considering frequency. 

Polarization has no pronounced effect on fhe image brightness from open water 
except at near range depression angles (Figure 9-10; A ) with usual over-land gain 
settings because returns of both polarizations appear black. However if the gain 
were set high enough to image signals from the water surface, W signals would be 
much brighter than HH or cross-poiorized signals. Marsh areas (m) dominated by 
the Black Spruce and wetland regions (w) dominated by grasses (Figures 1 1-12) 
exhibit strong variations in the degree of returned energy for different polarizations 
(HH or HV) at Ka-band frequencies. One area (Figures llA and IIB; Areas X) 
exhibits tonal reversals opposite to that observed for the wetland regions and mareh 
areas. On like (HH) polarized radar imagery (Figure llA; Areas X), fhe radar 
refurn Is comparable fo fhe surrounding deciduous forest but on cross (HV) polarized 
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FIGURE 9 . AN/APQ-97 radar Imagery showing polarization effects on open 
water In the extreme near range. 
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FIGURE 10. AN/APQ-97 radar imagery showing polarization effects on open 
water in the extreme near range. 
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AN/APQ-97 radar imagery depicting tonal 
and cross (HV) polarization of wetland (w) 


reversals for like (HH) 
and marsh (m) regions. 
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FIGURE 12 . 


AN/APQ-97 radar imagery depicting tonal reversals for like (HH) 
and for cross (HV) polarization of wetland (w) and marsh (m) regions. 
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imagery (Figure IIB; Areas X) these same locales exhibit high radar returns. The 
vegetation is principally sawgrass which belongs to the Sedge family, Cyperaceae, 
and is thought to be of the large and complex genus Carex . It is postulated that 
the cross 0^V) polarized radar signal undergoes volume scattering, thereby 
producing rfie anomalous high return. This strong depolarization occurs only for 
areas dominated by the Carex . Presumably this relates to the tendency of Carex 
to "droop" when standing in water, while the other species observed remain 
more erect. 

It has been noted by previous investigators [MacDonald and Waite, 1971] 
that like-polarized imagery provides maximum differentiation at K“band wave- 
lengths; Baflivala andUIaby [1975] observed that polarization had little effect 
on the moisture response for bare ground for low microwave frequencies between 
2 GHz and 8 GHz. However, wetland areas and marshy conditions were not 
considered in the latter investigation. Operating at frequencies lower than X~band 
for marsh delineation and for estuary and wetland surveys may or may not require 
polarization consideration . 

Thus far, the discussion has centered on free-standing water such as reser- 
voirs or ponds, marshes and wetland regions, and flood inundated areas which are 
not covered by forest canopies. Pools of water beneath vegetation and particularly 
under canopies of dense vegetation cannot be detected with either X“band or 
L-band radar for moderate and low depression angles [Drake et al ., 1974] . More- 
over, several studies at K~band frequencies [MacDonald, 1969; Wing, 1971; 

Moore, 1971] indicate that penetration does not occur and that the return Is from 
the top of the tree canopy. Barr and Miles [1970] observed higher radar returns 
from vegetation, especially low vegetation, peripheral to pools of water and con- 
cluded that it was the result of vegetation with a higher moisture content, which 
is consistent with controlled observations of vegetation backscatter [Ulaby, 1975]. 
Thus it may be possible to indirectly detect pools under dense vegetation. 

It is apparent that detection of water under a dense canopy is impossible at 
moderate and low depression angles; steep depression angles, where the angle of 
incidence approaches 0^, may provide necessary penetration capability depending 
upon the canopy characteristics. For example, by increasing the wavelength from 
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0.86 cm to about 20 cm (Ka"band to L"band)/ the attenuation is decreased and there" 
fore greater penetration of the overlying canopy may be possible. However, it must 
be remembered that as the angle of Incidence approaches zero, strong backscatter is 
recorded from smooth surfaces which not only can create the interpretation problems 
as discussed for Figure 7, but may also result in the similar return intensities originally 
observed at moderate and low depression angles* Thus, another trade"off appears to 
exist when considering depression angle. 

A final category for discussion concerns flowing water. Much information can 
be gathered on sub-surface topography by noting the surface flow characteristics. 

Radar may have a potential application if the system could be designed such that 
small scale flow patterns could be observed, A low frequency system would probably 
be unacceptable since it would tend to smooth the very features necessary for any 
analysis. An adequate resolution cell should be equivalent to about one”fourth of the 
width of the river. For major drainage systems, 50 meter resolution is not unrealistic; 
for smaller rivers a resolution of about 10 meters should provide adequate detail , For 
shoreline analysis including deltaic phenomena, the latter resolution requirement should 
provide the maximum amount of detail required for most interpretations. 

4.0 CONCLUSIONS 


The optimum system parameters are difficult to establish at this time since 
little research has been undertaken to provide adequate data on sensor responses to 
land"water boundaries under differing physical conditions. 

Incidence angle (depression angle) considerations are an important parameter 
from the standpoint of signal return from open water. To insure accurate shoreline 
delineation and eliminate the problem illustrated in Figure 7, a limitation in the 
depression angle should be established at 75°. This limitation may or may not provide 
the necessary penetration capability needed for pool detection under vegetation 
canopies. A major control parameter will be vegetation density. 

An average revisit time of not more than 6 days is suggested as a compromise 
for lakes, ponds, and estuaries (size of feature and type of data requested obviously 
control periodicity of coverage needed); wetlands require less frequent coverage, 
from 30 days to 1 year depending upon their physical size. For flood plain delinea- 
tion, channel characlerisHcs, etc., periodic coverage each year would suffice, while 
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flood damage assessmenf requires coverage upon demand (periodicity of coverage 
obtained from tables of measurable parameters from a report entitled “Satellite 
Data Collection User Requirements: Report of the NASA Ad Hoc Committee on 
Hydrology"). Thus as illustrated above, revisit time can fluctuate depending upon 
the feature for which data is recorded. To insure adequate periodic coverage for 
all surface water phenomena, a revisit cycle of 6 days is suggested as a compromise, 
but a system should be capable of "pointing" to flooded areas on demand. 

Polarization is an important parameter for water surfaces imaged at steep 
depression angles, the extreme near range. Cross"polarized imagery (HV) however, 
subdues the radar response to marshy areas. The least complicated solution appears 
to be the utilization of like-polarized ^H) imagery at depression angles shallower 
than those generating the problem response in the near range. However, the use of 
dual polarization (HH and HV) allows for better identification of vegetational 
differences within wetland regions. 

Based upon existing radar imagery viewed by this investigator, a proposed 
system resolution of approximately 12 meters should provide sufficient information for 
both regional analyses and provide for local detail, although finer resolution would 
produce greater detail , Additional data should be acquired however to quantify this 
requirement; one such investigation was conducted by Moore and Dellwig [1976] for 
a variety of target types . 

Frequency will control the radar responses observed for any given target, 
past investigations have dealt primarily with X“band and K-band imaging systems. 
Based on the limited research data base existing for water body detection, L“band 
does not appear to possess advantages over higher frequency systems and in many 
cases, as previously discussed, resluted in serious limitations. Unfortunately, data 
is not available for S“band or C”band imagery, thus it is impossible to determine its 
feasibility for monitoring the features discussed above. 

Based on available research, the following recommendations are proposed: 


frequency: 

X-band 

polarization: 

HH 

depression angle 
limitation: 

less than 75 ' 

resolution: 

12 meters 

revisit time: 

6 days 
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In riie opinion of the author, a system operating at 4 GHz (proposed system for 
soil moisture) should provide the necessary data without the disadvantages observed 
at L“band frequencies. This contention cannot, however, be substantiated with sup- 
portive data. Research is limited to a very few investigations at high frequencies and 
even fewer which utilize the lower frequencies. The obvious solution to the dilemma is 
research aimed at the whole surface water problem; research that will provide data for 
frequencies between L-band and K-band such that meaningful comparisons can be 
realized. A working ground-based research system exists whereby experimental data 
is retrievable for frequencies between 1 GHz and 18 GHz and over all the potentially 
useful incidence angles and polarizations; a major obstacle delaying such an investigation 
is project funding . 
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ABSTRACT 


The potenHal for monitoring ice by imaging radar on the Great takes has 
been demonstrated, but other applications of radar to monitoring of ice, for instance 
in small lakes and ponds of hydrologic significance, remains to be explored. Although 
ic6*“covered rivers appear in numerous radar images, no systematic study of the river- 
ice monitoring problem has yet been made. The only frequencies studied in the Great 
Lakes experiments have been in the L" end X“bands, with most measurements in the 
latter band. No systematic study has been made of the effect of angle of incidence, 
polarization, or resolution. Clearly considerable experimental work is needed in this 
area. A theoretical model has been developed for sea ice, but only theory-based 
Inferences have been used to describe radar-return from lake Ice. 

The formation of lake ice is reviewed with its implications for, and some 
measured data for, the electrical properties. At the lower frequencies penetration 
in ice should be great enough so that most of the return probably comes from the bottom 
of ice in the Great Lakes, Suggestions for Ice-monitoring radar system parameters are 
made, but they are based largely on conjecture from the physical and dielectric char- 
acteristics of the ice. More important at this stage are the recommendations for a 
research program. 



Radar Monitoring of Lake/River Ice “ 
State of the Art 
by 

S. K. Parashar 


1.0 INTRODUCTION 

The current state of the art in monitoring lake and river ice by the use of a 
radar is presented here on the basis of a thorough review of the existing literature. 

The present Interest in the use of radar systems to monitor lake and river ice was 
primarily developed because of a real need to extend the Great Lakes shipping season 
into the winter months. The feasibility of achieving this depends to a large extent on 
making improvements in the ice information gathering techniques. It is important that 
quick, accurate, and comprehensive information about the position, extent, and rela- 
tive thickness of ice cover be made available to shippers on a timely basis so that ship- 
ping routes for navigation may be optimized. This requires repeated, and in some coses 
daily, reconnaissance of the ice cover. In view of the tremendous areal extent of the 
ice cover, the repeated and timely surveillance can be only provided by means of auto- 
mated remote sensing techniques. It is in this regard that the all-weather, day/night 
operational capability and broad aerial coverage provided by a Side-Looking Airborne 
Radar (SLAR) is seen to be of great potential utility in the ice mapping. 

Besides o comprehensive review of the past radar measurements of lake and river 
ice, a section on the formation of lake and river ice is included and the physical and 
electrical properties of ice which help determine the radar return are given. This will 
not only help in understanding the nature of radar return from ice but also help in 
explaining the experimental measurements. In the end an attempt is made to specify 
optimum parameters for an operational ice surveillance system on the basis of past 
experimental results. Recommendations are made on the need to conduct future experi- 
ments and on the design and requirements of these experiments, 

2.0 FORMATION OF LAKE/RIVER ICE 

The formation of ice depends primarily on the surface salinity, the vertical 
distribution of salinity, and the depth of water. The amount of salinity present deter- 
mines the temperature of maximum density. In water with a salinity of less than 24.7 °/oo 


the temperature of maximum density lies above the freezing temperature. The freezing 
temperature of water decreases less rapidly than the temperature of the maximum den" 
sity, with increasing salinity. The temperature of maximum density in pure water is 
+3,98° C. Both the temperatures, of freezing and maximum density, are the same at 
-1 ,33° C at a salinity of about 24.7 '^ 00 . Cooling of a layer of water with a salinity 
of less than 24»7 ^^oo causes an increase in the density of water at the surface resulting 
in free convection. With continued cooling this convection continues until the complete 
convective column reaches the temperature of maximum density for that particular 
salinity (14), The convection stops at this time. Further surface cooling, without 
mechanical mixing, creates a stratification in a very thin stable surface layer of low 
density where ice formation rapidly begins when the temperature of the freezing point 
is reached. Lake or river water is of considerable less salinity than 24,7 *^ 00 . 

The formation of ice at the surface begins earlier over shallow depths than over 
deeper water bodies under similar conditions (5), In the fall lake ice normally forms 
first, then river ice, and finally sea ice. Both lake ond river ice are essentially fresh 
water structures whereas sea ice is quite salty. The amount of salt present in ice depends 
on the salinity of the water from which it is formed and is always less than the original 
water salinity. Ice formation in the rivers is affected by the complete vertical mixing 
of heat as the result of turbulence. 

The salinity of the water does not affect the formation of the initial ice cover 
once the water starts to freeze. The amount of turbulent mixing of the water in the 
freezing layer is the determining factor in the formation of ice cover. Thus, factors 
such as wind, current, and the Intensity of cooling are important in influencing the ice 
formation. Under calm conditions, the initio! ice growth starts after water has been 
supercooled slightly. The first crystals to form are minute spheres of pure ice and as 
these grow they rapidly change their shape to circular disks (14). In fresh water these 
disks attain a maximum diameter of 2 to 3 mm which depends on supercooling. The disk 
form changes to dendritic hexagonal stars at some critical diameter of about 2 to 3 mm 
in fresh water. This critical diameter appears to decrease with increasing salinity. 

The stars grow rapidly across the surface of calm lake water until they overlap and 
freeze together. 

Surface needles are common during the initial freezing of fresh water and are 
formed when a disk becomes inclined at an angle to the water surface. Subsequent 
growth is in the form of long thin needle~Iike crystals. These crystals vary considerably 
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in size and needles up to 4 m in length and vertical crystals up to 0.7 m in diameter are 
observed (14). Thus, during calm conditions the initial lake ice skim is composed of an 
open polygonal network of surface needles inclined at some angle from the vertical and 
in between areas are occupied by the approximately vertical starHike crystals. 

Some turbulence during initial ice formation introduces more nuclei into the 
area of active freezing so that the initial crystals are mixed throughout a depth of up 
to several meters. Extensive discordal growth is favored and in rivers these discords 
and needle"like ice fragments are referred to as frazil ice. Irregular structures are 
formed when these frazil crystals freeze together. Snow falling on the water surface, 
during calm conditions, just prior to and during initial freezing conditions also helps in 
the formation of similar crystal aggregates. Ice known as pancake and slush is formed 
because of packing of ice crystals due to wave motion and the general motion of the 
aggregates against one another. The thickness of initial ice cover when a slush layer 
congeols is usually several centimeters. Sheet ice will develop beneath the pancakes 
and slush ice once a composite ice sheet is formed between the pancakes. It is difficult 
sometimes to distinguish slush Ice from infiltrated snow ice which is formed because of 
abundant snowfai I , 

After the initial continuous skim of ice is formed, a transition layer is produced 
in which a change towards preferred growth orientation of ice crystals is seen. The 
favored orientation can range from horizontal to vertical. The ice below the transition 
layer has all the characteristics associated with a ‘'columnar" structure. This columnar 
zone has a strong crystal orientation and a gradual increase in average grain size as the 
distance from the surface or the cold source increases. Several characteristics and 
features such as air bubble layering can be seen in the vertical sections of columnar 
zone. The shape of the air bubbles is usually tubular and elongated parallel to the 
growth direction. In lake ice, a pronounced intra“crystalline substructure, which 
usuallv forms parallel to the growth direction and has a similar appearance to so called 
striatlon boundaries In metals, can be seen. 

After the formation of initial ice skim, both lake and river ice usually grow 
with a planar solid/liquld interface. It is pointed out by Weeks and Assur (14, that 
the voriations in the properties of lake and river ice with changes in temperature should 
be caused primarily by variations in the properties of pure tee and not by any changes 
in the relative amount of solid, liquid, or gas present in the Ice sample. 
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2,1 Physical Properties 




Both lake and river ice consist of pure ice crystal matrix in which air bubbles 
are included. The ice may also contain certain forms of solid impurities and liquid 
inclusions depending on the state and purity of original water and the conditions pre“ 
valent at the time of ice formation. The properties of lake and river ice will be 
primarily those exhibited by pure ice. It is possible to classify natural Iy"occuring 
fresh water ice into two major categories on the basis of its formation and growth: 

(1) ice which grows parallel to the direction of heat flow is called secondary ice; 

(2) ice which forms on top of the secondary ice, for example, as the result of the flood 
ing of the snow layer, !s termed superimposed ice. The secondary ice can be further 
classified into several ice types on the basis of age and thickness and the v/ay it is 
formed. From the point of view of remote sensing, frazil Ice and columnar ice are 
probably most interesting (13). As pointed out earlier, frazil ice Is formed from an 
agglomeration of ice particles which have been formed in supercooled turbulent water 
and can be several centimeters thick. Columnar ice is formed in calm water or under 
an already existing ice cover, such as frazil ice. 

It is the physical properties of ice such as amount and nature of impurities, air 
bubbles, cracks, density, and temperature which determine the electrical properties. 
The radar return from lake and river Ice depends on the surface roughness at the wave" 
length scale and the subsurface structure. 


2.2 Electrical Properties 

The radar return from lake and river ice, as from other surfaces, is determined 
by the electrical and physical properties of ice. The electrical properties depend on 
physical properties. The electrical properties of ice of Interest ore Its dielectric be" 
haviour such as dielectric constant and conductivity as a function of ^ o uency and 
temperature. 

The electrical properties of fresh water or pure ice have been investigated by 
many workers and are quite well established. An overall review of literature on elec" 
trical properties of fresh water ice presented by Evans (4). As liquid water changes Its 
physical state into ice, its dielectric properties change sigr\iticant!y . The most Inter" 
esting properties of pure ice are its high static dielectric constant (about 100) and Its 
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long reiaxafion time (about 10 For frequencies much greater than 1 MHz, the 
dielectric constant drops to about 3. Most workers have investigated the laige diei" 
ectric dispersion in pure ice. In general, measurements of contemporary workers have 
a good correlation with each other and with the Debye equation. The Debye equation 
for dielectrics is based on polar molecules having a single relaxation time and this 
relaxation fomi is fairly common in solids and liquids. 

The dielectric behaviour of ice and woter as a function of frequency is given 
in Figure 1 . Both ice and water traverse in first approximation a simple Debye reIoxa~ 
tion spectrum. In the dispersion region the dielectric loss goes through a maximum and 
the dielectric constant falls off. The separation in the relaxation spectra of ice and 
water is approximately six decades of frequency (6) . The dielectric loss of water and 
ice changes drastically by the presence of salts whereas the dielectric constant does 
not alter as much. Salts add free charge carriers and thus increase the conductivity. 
Thus, the electrical properties of sea ice which contains salt are significantly different 
from those of pure ice, and lake or river ice is more like pure ice. 

Both the relaxation time and the static dielectric constant of ice are temperature 
dependent. Between a frequency of 1 MHz and the far Infra-red region theie \s no 
absorption band in the spectrum of ice. The relative permittivity (dielectric constant) 
from the measurements made by different workers and compiled by Evans (4), is shown 
In Figure 2, The code letters identify the individual authorities along with the tempera- 
ture, and other remarks about the measurements. The differences in the measurements 
can be attributed to differences in the temperatures and densities of the tee samples. 

A scatter plot of the experimental results of dielectric constant obtained by 
various samples of fresh ice as a function of temperature at 10 GHz is shown in Figure 
3. Statistically the value of relative dielectric constant was found to be 3,14 + 1,4%, 
which concurs with the measurements made by other workers (13), For fresh (tap) water 
Ice, the measurements in the 26.4 to 40 GHz frequency range indicate a value of rela- 
tive dielectric constant of 2,92 + 2%, 

The variation of loss tangent with frequency for ice as compiled by Evans (4) 
is presented In Figure 4. The quantify plotted vertically is log^Q (f tan 5} where f is 
the frequency In MHz and tan 5 rs the loss tangent or € ^ / ^ (imaginary part of 

complex permittivity divided by the real part). 

The loss tangent of fresh ice at 10 GHz as obtained by Vant et ol (13) was of 
the order of 20 X 10 {0° C to -35° C). Lamb (13j gives a value of 12X10"^ 
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(0® C), Cummings (13) 27 X 10 ^ (0° C) to 6,5 X 10 ^ (-18°C), and von Hippel (13j 

7.0 X 10“"^(-12° C), 

The presence of impurities in small quantities may have some effect on the 
dielectric properties. This effect may be significant depending on the amount and the 
nature of impurities. 

3.0 RADAR RETURN FROM ICE 


The radar scattering from lake and river ice is a result of complex interaction 
of the electromagnetic field with the ice profile. The average power returned by each 
boundary (air-ice and ice**water boundary) is partly determined by the complex permit- 
tivity and the surface roughness. In addition subsurface structure may also have some 
effect on the radar return. A perfectly smooth surface gives return only in the specular 
direction where angle of incidence is equal to angle of reflection. With a very rough 
surface, the scatter is almost uniform in all the directions including the source. For 
the slightly rough surface most of the incident energy is scattered in the direction near 
the specular direction. The surface roughness Is defined on the incident wavelength 
scale. For o surface to be defined as smooth, the surface variations should be much 
smaller than one wavelength, perhaps by an order of the magnitude. Thus, a surface 
which is smooth at one frequency may not be smooth at a higher frequency. The same 
surface will appear rougher as the wavelength of the incident energy is decreased. 

In addition to surface roughness, the electrical properties of the medium also 
play an importont part. The complex permittivity (the real part is the dielectric con- 
stant anJ the imaginary part can be expressed as a loss tangent) not only determines the 
amount of energy which is scattered or reflected but also determines the degree of 
penetration of incident electromagnetic energy. The magnitude of complex permittivity 
determines the amount of energy scattered or reflected and the loss tangent determines 
the attenuation in the medium and thus the skin depth and the degree of penetration. 

For higher values of the loss tangent, the attenuation is greater and hence the skin 
depth smaller. By selecting a value of 3,15 for dielectric constant and 10 X 10 
for loss tangent, the skin depth for ice can vary linearly as function of frequency from 
about 100 cm at 1 GHz and 10 cm at 10 GHz to 1 cm at 20 GHz, 

Thus, the choice of the frequency determines the depth of penetration and there' 
by influences the contributions made by the ice-water interface and the effect of sub- 
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surface structure. The contribution to radar return from surface roughness is also 
dependent on frequency. 


4.0 STATE OF THE ART 


The present state of the art in the radar measurement of lake and river ice 
con essentially be obtained by conducting a thorough review of the existing literature, 
and this is the intent and purpose of this section. 

Most of the past research efforts in the applications of radar in monitoring ice 
have been devoted to sea ice. The investigations into the ability of radar to monitor 
lake and river ice have been rather limited until quite recently. This is partly because 
of the fact that it is felt that the knowledge about sea ice is of greater practical utility 
and global Interest than the knowledge about lake ice which is of regional valuf , It 
is only recently that interest has been developed in investigating the potential of 
radar in monitoring lake and river ice, 

r During February and March, 1971, on attempt was made by United States 
Coast Guard to acquire some SUAR {SIde"Looking Airborne Radar) imagery in the 
Mackinac"Sault Ste, Marie region of the Great Lakes, The radar imagery was obtained 
with an AN/DPD-2 (Modified) Philco-Ford SLAR operating at 16.5 GHz:. In addition 
to SLAR imagery, vertical aerial photographs in the 9X9 inch format were also obtained 
of certain areas. This was done so that radar imagery could be correlated with the 
aerial photographs. It was shown then by Photographic Interpretation Corporation in 
a report prepared for U.S, Coast Guard (11) on the analysis of the data that radar 
imagery is valuable for the interpretation of lake ice. It was possible to detect, delin- 
eate, and describe several lake”ice types through a detailed, systematic study of the 
available imagery. This was accomplished even though there was lack of information 
available at that time concerning the various basic categories of lake ice types and 
their parameters (e.g., thickness, age, and strength). At Ku“band no apparent detri- 
mental masking effect was found by snow cover of lake ice. Some instances of prob- 
able penetration of snow cover by radar signals were observed and described. No snow 
features identifiable as such were detected. Some distinct patterns on the radar imagery 
were attributed to snow, ice ridges, and mounds on the lake ice surface. It was pos- 
sible but difficult to describe lake ice patterns and determine relative thickness. 

The same data were examined by Raytheon Company in a report prepared for 
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U.S. Coast Guard (12), It was pointed out that identification of the various ice 
types was accomplished to the maximum extent possible. It was possible in certain 
instances to identify new ice types with the aid of complementary photography. 

Slush, frazil, and grease ice sometimes were not differentiated on lake ice imagery. 

It was not possible generally to separate young ice into dark gray or gray"whfte types; 
these had to be considered one unit. It was only possible to separate dark nllas and 
light nilas with the aid of aerial photographs. Not enough ice rind was found on the 
SLAR imagery to arrive at definite conclusions for discrimination, A glossary of ice 
terms generally used is provided by Dunbar (3), 

Winter ice, considered to be of "medium" thickness, was interpreted with 
relative ease. It was pointed out that even qualitative thickness could not generally 
be determined with SLAR to any degree of confidence. In certain cases, winter ice 
appeared to be thicker, and if the surrounding features and collateral data supported 
the assumption, such ice was termed "thick" winter ice. Sometimes, clues to relative 
ice thickness were provided by the crack system. Angular cracks implied thinner ice. 

This technique is not a reliable one, ^ 

It was pointed out that the masking effects of snow cover on lake ice tended to 
complicate the identification of ice types. As confirmed by aerial photographs, almost 
all of the ice having a continuous snow cover was winter ice. Still other variables are 
introduced by thin and/or discontinuous snow cover which has to be considered in the 
identification. If was recommended thot the effects of such lesser amounts of snow 
cover be analyzed further under controlled conditions so that these influences on the 
SLAR imagery may be understood and isolated, 

A series of X"band SLAR images were obtained by NASA Lewis Research Center 
(7) to show the development and disintegration of the entire ice cover on Lake Erie during 
the winter of 1972-73. The accurate correlation of radar responses with ice conditions 
was established through simultaneous ground truth observations and ERTS-1 photography. 
Motorola AN/APS-94C real aperture SLAR system operating at 9.2 GHz frequency 
with horizontal transmit and receive polarization was used to acquire imagery. It was 
possible to identify ice types such as brash, pancake, and related forms because of their 
brightest return. It was pointed out that these ice types give brightest return simply 
due to the large vertical cross-section presented by their edges, Thusj, surface rough- 
ness of ice was considered to play a dominant role in the radar return. Because of the 
penetration of the wave the roughness of the ice“Wdter interface may also contribute 
strongly to the radar return. Only in the case of fast ice were there any indications of 
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volume scattering. It was not possible to discriminate unfractured ice from open 
water because {it was felt) smooth clear ice lacks sufficient defects to backscatter 
enough radiation to be detected. Most of the returns in ihis stud/ were near 
grazing incidence, so the effect of steeper angles cannot be evaluated. 

During the 1973-74 winter, measurements were made with the 2-frequency 
(X-band and L-band) synthetic-aperture system of Environmental Research Institute 
of Michigan (2 ), Two sites near the entrance to Great Lakes were studied: White- 
fish Bay on Lake Superior and the Straits of Mackinac between Lakes Michigan and 
Huron, Both HH (horizontal transmit, horizontal receive) and HV (horizontal trans- 
mit, vertical receive) images were produced in each band. The systems have resolu" 
tions much finer than for other reported measurements. 

Smooth black ice with embedded brash ice gave weak returns and rough 
brash ice gave strong HH returns at both frequencies. In one relatively smooth 
area moderate return occurred only on the L~band (HH, HV) images. One area 
of poorly-developed ice foot showed only on X-band (HH, HV) images. In most 
oases pressure ridges showed up on all images, but one could only be detected on 
the X-HV image. 

Interpretations suggested that much of the scatter was from the lower surface 
of the ice, particularly at L-band. 

As is clear from above, the radar measurements from lake ice are rather 
limited and there are no measurements available for river ice. Furthermore, the 
past results are inconsistent. The only consistency which can be seen from differ- 
ent studies is that radar may very well be a valuable too! in discriminating and 
Identifying lake ice. It is not possible to measure ice thickness directly, but it 
can be inferred by dividing ice into types on the basis of thickness. There is a 
lack of qualitative and quantitative data available which can help in relating ice 
thickness and/or type to radar backscatter. More research should be conducted 
with a view not only to explore the potential and prove the capability of SLAR in 
monitoring lake ice but also to provide enough reliable information for the design 
of future operational systems, 

5.0 RECOMMENDATIONS 


To be able to design an optimum operational radar system for ice monitoring, 
the effect of frequency, polarization, and angle in discriminating ice types has to 
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be clearly u. derstood. This understanding of the nature of radar return can be 
obtained in part by collecting qualitative data corresponding to different operating 
parameters and then establishing the optimum ones. This, however, requires much 
data gathering and the time required to accomplish this and the cost may be pro“ 
hibitive. 

Another approach is to formulate an analytical model which will help In 
computing the radar return corresponding to different operating parameters and ice 
types and conditions. For this case, the qualitative data with excellent ‘'ground- 
truth" need only be acquired for selected parameters so that the validity of the 
theoretical model can be established. On the basis of this model optimum paramet" 
ers can then be specified. Such a theoretical model has been formulated for sea 
ice (10) and it shows reasonably good agreement with the experimental results. The 
good point in the favor of a theoretical model is not only that it helps in explaining 
and understanding the nature of radar return but also it separates out the effect of 
surface roughness from other effects. Thus, the model helps in establishing the 
relative contributions from surface roughness and the dielectric properties. 

In addition to establishing the frequency, polarization, ond angle, the 
other parameter which needs to be established is the size of the resolution cell. It 
may be possible to discriminate ice types better with a moderate resolution cell 
size than with fine or coarse size. It is not suggested here that only a single fre- 
quency, polarization, and range of angles is suitable for ice discrimination. It 
may very well be that a combination of frequencies, polarization, and ranges of 
angles provide better information in identifying ice types. 

The other parameters which need to be specified are the sensitivity and the 
dynamic range of the radar return. This can be achieved by acquiring radar scatter- 
ometer measurements over varying ice conditions. Radar scatterometers by design 
are calibrated systems and they measure the radar backscatfer coefficient or radar 
scattering cross“section per unit area. Radar scatterometers, by measuring voriation 
of the radar scattering coefficient, o-°, with angle, frequency, and polarization, 
permit more detailed observation of radar scattering behavior than radar Images (9), 
Almost all the radar imagers in operation today are uncalibrated systems operating 
largely near grazing incidence. The gray tone or density on a radar image is pro- 
portional to the power received which in turn is directly proportional to the radar 
scattering coefficient. Thus, by establishing the sensitivity and the dynamic range 
of the radar backscotter coefficient the corresponding parameters for a radar Imager 
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can be specified. No radar scafi'erometer data is avilable for lake and river ice. 

Though it is a difficult and a dangerous task to specify parameters for future 
operational systems on the basis of meagre data available/ yet an attempt is made 
below to infer these parameters. The specifications recommended here are by no 
means optimum or the best or the only ones possible, 

1) Frequency 

Past measurements have only been made at two frequencies: X“band and Ku” 
band plus a very few of L"band, It does not appear that one high frequency is better 
than the other, though more ice types were identified on Ku”band imagery. In 
certain instances the snow masking effects were detrimental to ice identification but 
in general Ku“band signals do appear to penetrate dry snow. 

The two parameters of interest in ice monitoring are the relative ice thickness 
and the surface roughness. Both of these parameters help determine the mechanical 
properties such as strength and also help in establishing the drift patterns and climato" 
logical conditions from the corresponding mathematical theoretical models. The 
knowledge of ice thickness and also the presence of open water helps in determining 
the optimum navigational route. 

The choice of frequency or frequencies should then be mode by bearing in 
mind the above mentioned parameters of interest namely. Ice thickness and rough- 
ness. A signal at a lower frequency will generally be able to penetrate more as com- 
pared to a higher frequency. Thus, radar return at a lower frequency will be 
relatively more influenced by the ice-water boundary. Moreover, the contribution 
mode by the roughness is also dependent on the frequency. The same surface will 
appear rougher with an increase in frequency. 

The frequency or frequencies should also be selected on the basis of achieve- 
able resolution cell size with a reasonable size antenna. 

By keeping these factors in mind, a choice of two frequencies, one at L"band 
around 2 GHz and the other at Ku-band around 16 GHz seems to be appropriate. 

If only one frequency is required, X"band around 10 GHz seems to be more feasible. 
The signals at the L-band will be able to penetrate more and thus the radar return 
will be influenced to a greater degree from the ice"water boundary. This may help 
in differentiating thicker types of ice. Also L-band signals will be effected more 
by the large scale roughness and may very well provide qualitative and quantitative 
information about the roughness such as pressure ridges. Ku“band signals will be 
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able io penetroi-e less and will be more influenced by the small scale roughness. 

This may help in differentiating thinner, very rough types of ice such as brash. 

The choice of X"band is made as single frequency simply because conceivably 
it will have the best of both. It will be able to penetrate a little and also be in- 
fluenced by medium scale roughness. To achieve moderate or fine resolution at these 
frequencies, especially at LHsand, a synthetic aperture antenna is required even 
for aircraft systems. SA.R is required for spacecraft systems at any frequency, 

2) Polarization 

There does not appear any difference between two like polarizations (HH - 
horizontal transmit, horizontal receive; W - vertical transmit, vertical receive) in 
their ability to distinguish ice types, but cross-polarized signals may be effective in 
discriminating ice types as is evident in the case of sea ice {10). Thus, one like 
polarization (HH or W) and one cross-polarization (HV or VH) is recommended, 
but the latter is based largely on conjecture. 

3) Incidence Angles 

The selection of incidence angles should be based on their ability to distinguish 
ice types. No such Information is available from the past measurements. Since one 
of the aims in Ice monitoring is the location of navigational open water, the choice 
of angles may very well be based on this criterion. Open water can generally be 
detected by operating well away from the vertical because as a comparitively smooth 
surface it does not give any return at near grazing angles and thus appears black. 
Hence a range of angles from about 30 to 75 degrees is proposed. Since some of 
the smooth ice categories also give low return, these angles may be less suitable 
than they are for sea tee. 

4) Resolution Cell Size 

This selection should be based on two criteria. One Is the minimum resolution 
cell size needed to identify on ice type and the minimum size floe which needed to 
be identified. The other is the minimum size of open water lead which is adequate 
enough for the passage of a ship. For lack of better information available, a mod- 
erate size of 20 to 30 meters is proposed. Note that this is finer than that of either 
DPD-2 or APS-94C, 

5) Dynomic Range 

No data are available on which selection can be based. On the basis of sea 
ice data (10), a dynamic range of 40 dB at above mentioned operating angles should 
be more than adequate. A dynamic range of less than this may actually be required. 
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6) Sensitivity 

Again there is a lack of available information; in fact there is none! On 
the basis of sea ice data (10), it should be possible to detect a minimum radar scatter 
ing coefficient, of about “20 dB for like polarization and "50 dB for cross- 
polarization* The system should be sensitive enough to record and distinguish var- 
iations of 0.2 dB and these variations should be reflected in gray tones on the radar 
imagery. Some evidence indicates that sensitivity better than for sea ice would be 
useful . 

7) Frequency of Coverage 

During winter time a coverage every day may be needed to find optimum 
shipping routes on the lakes. Coverage for rivers may also be needed daily at 
times of ice-jum build up, but weekly coverage may be sufficient at other times. 

6.0 CONCLUSIONS 


It Is evident from the above that SLAR is potentially a valuable tool in mon- 
itoring lake and river ice. Before an optimum system of general utility can be 
designed more research, both experimental and theoretical needs to be conducted. 

A need exists to collect more qualitative and quantitative dote with a view of 
establishing optimum operating parameters and understanding the nature of radar 
return from ice. These experiments should be carefully designed and properly organ- 
ized so that maximum information from them can be derived. The need to collect 
good '‘ground truth" data in terms of ice thickness, surface roughness, temperature, 
and electrical parameters should be borne in mind. The experiments should be 
conducted in different parts of the season so that the effect of temperature, growth, 
and melt can be properly evaluated. The measurements corresponding to different 
frequencies, polarizations, range of angles, and resolution cell size should be made 
for both lake and river ice. A "ground-based” radar spectrometer offers more 
potential for these measurements than any airborne system. 

A theoretical model of wave scattering from fresh-water ice needs to be 
formulated. To test this model properly, measurements under semi -control led 
conditions have to be made at selected parameters. To properly correlate these 
measurements with ice types a good system to collect the necessary "ground truth" 
information must be set up. These measurements can be made by a ground based 
scatterometer operating at a number of frequencies. Later airborne scatterometer 
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measurements may be helpful to establish dynamic range and sensitivity. 

The effect of snow on radar return needs to be evaluated more thoroughly. 
Furthermore, lake and river ice should be divided into standard categories as has 
been done in the case of sea ice* Relative thickness should be associated with 
each ice type. 

Fine "resolution images should be produced to assist in establishing required 
resolutions. These may be degraded by different amounts, and the degraded images 
evaluated by interpreters to set critical levels for different purposes. 
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Figure 1 . The dielectric behavior of ice and water as a function of 
frequency {Hoekstra and Capillino (6)), 
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Figure 2. Relative permittivity of ice (ordinates) versus logarithm of 
radio-frequency (abscissae) (Evans (4)). 

L: Lamb (1946) and Lamb and Turney (1949) - 5° C at low frequencrs, 0° to 
-190 C at high frequencies: distilled water. 

C: Cummings (1952) "18 C. Distilled water and melted snow. 

A: Auty and Cole (1952) - 10 C. Conductivity water: ice free from stress. 

V: Von Hippe! (1954) -^12 C.^Conductivity water: ice not annealed. 2 

Y: Yoshlno (1961) “■ 18 to -36 C. Antarctic ice^ not annealed, density 0.91 g/cm. 
W; Westphal (priva^ communication) “5 to “60 C., annealed Greenland ice, 
density 0.90 g/cm. 
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Figure 1 . The dielectric behavior of ice and water as a function of 
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Figure 2. Relative permittivity of ice (ordinates) versus logarithm of 
radi ©"frequency (abscissae) (Evans (4)). 

L: Lamb (1946) ond Lamb and Turney (1949) “ 5° C at low frequencis, 0*^ to 
”190“ C at high frequencies: distilled water. 

C: Cummings (1952) ”18 C, Distilled water and melted snow. 

A: Auty and Cole (1952) - 10 C, Conductivity water: ice free from stress. 

V: Von Hippe! (1954) “12 C. Conductivity water; ice not annealed. o 

Y: Yoshino (1961) - 18° to ”36 C. Antan:tic ice^ not annealed, density 0,91 g/cm, 
W; Westphal (private communication) “5 to "60 C,, annealed Greenland ice, 
density 0.90 g/cm. 
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Figure 3« Results of free'space rneosur'^ments of dielectric 
constant of pure ice (Vont, et al., (13)). 



Figute 4. Lpito tangent of led versus rqdio^&equency. The quantity plotted 
veiticoliy is log|Q (f ton 6 ) where f is the frequency in Mc,/sec. On the 
high-frequency ' tail of o reioxation spectrum this quontity is constant: it 
has the further useful property that the attenuation of a radio wave (measured 
in dB./m.) passing through trie medium is directly proportioned to f tan 5 . 
Temperatures are marked in C. (Evans, (4))* (key on folfovy'ing page). 


L: Lamb (1946) and Lamb and Turney (1949), DtsHlIed water, ice not annealed^ 

C: Gumming (1952), Distilled water, tap water, and melted snow (no observable 
difference) , 

A; Auty and Cole (1952). Conductivity water, ice free from stress. Limiting 
values plotted arbitrarily at 1,000 times the relaxation frequency, 

V: Von Hippel (1954), Conductivity water, ice not annealed, 

Y; YoshinG'(1961) , Antarctic ice core samples, not annealed, density 0.91 g./cm. 

W: Westphal (private communication), Greenland ice, annealed, density 0,90 g./cm. 

Approximate temperature coefficients below “ 10® C. 

1 Mc./sec. 0,05 per ®C. in log tan 5 = 12% per ® C. In tan 6 (from Auty and 
Cole), 

iOO Mc./sec, 0.025 per ®G, in log tan 6 = 6% per ®C, in tan (S (from Westphal) 
Mc./sec. 0.01 per ®C. in log tan $ ~ 2.5% per ®C. in tan 6 (from Lamb) 
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